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ABSTRACT 

Recently the series for two RG functions (corresponding to the anomalous di- 
mensions of the fields 4> an d 4> 2 ) of the 3D cf> field theory have been extended to 
next order (seven loops) by Murray and Nickel. We examine here the influence 
of these additional terms on the estimates of critical exponents of the TV-vector 
model, using some new ideas in the context of the Borel summation techniques. 
The estimates have slightly changed, but remain within errors of the previous eval- 
uation. Exponents like r\ (related to the field anomalous dimension), which were 
poorly determined in the previous evaluation of Le Guillou-Zinn- Justin, have seen 
their apparent errors significantly decrease. More importantly, perhaps, summa- 
tion errors are better determined. 

The change in exponents affects the recently determined ratios of amplitudes 
and we report the corresponding new values. 

Finally, because an error has been discovered in the last order of the published 
e = 4 — d expansions (order e 5 ), we have also reanalyzed the determination of 
exponents from the £-expansion. 

The conclusion is that the general agreement between e-expansion and 3D series 
has improved with respect to Le Guillou-Zinn- Justin. 
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1 Introduction and summary of results 

Recently the perturbative expansions of the anomalous dimensions of the fields 
(ft and (ft 2 for the O(N) symmetric (</> 2 )^ =3 field theory have been extended to 



next order (seven loops) in the case iV = 0, • • ■ 3 by Murray and Nickel [ 4"6|j . This 
rather impressive result has led us to reexamine the determinations of the critical 
exponents for N = (polymers), N = 1 (Ising like systems), N = 2 (superfluid 
Helium) and N = 3 (real ferromagnets). For completeness we have added results 
(at six loops) for N = 4 which correspond to the Higgs sector of the Standard 
Model at finite temperature. A limitation of the present work is that the series 
for the RG /^-functions have not been extended (they remain at six loops) and 
for several exponents this now is the main source of error. 

Critical exponents have also been calculated in the form of e = 4 — d expan- 
sions, up to five loops jnj. Recently a slight error in the previously published 



series has been corrected ||12|1 , and this has motivated us to also reexamine the 
corresponding estimates (again adding N = 4 results). 

For the reader who is not interested in details we summarize our main results 
for N = 0, • • • 3 in Table 1 (d = 3) and in Table 2 (e-expansion) while N = 4 
results for both methods can be found in Table 3 . We have chosen central values 
which satisfy all scaling relations, but the apparent errors for 7, /3, r\ in general 
have been determined independently. For the d = 3 IR fixed point value g* we 
give results both in the usual field theory normalization (Eqs. (2.3)) and in the 
normalization used by Nickel , 

N + 8 

which is such that the fixed point value is close to 1. 

Note that in Table 1 in addition to the plain e-expansion results (denoted as 
"free") we report some additional results denoted as "be" (i.e. with boundary 
condition) that try to incorporate the knowledge of the exact d = 2 values by 
summing the series (f(s) — f(2))/ (2 — e), where f(e) is an exponent with known 
2D value. In the case of the exponent v for N = the d = 1 value is also known. 
We have checked that incorporating this additional piece of information has no 
significant impact on the final result. 

For N > 2 the analysis of the series with boundary conditions is quite difficult. 
Therefore we present here only central values, but no error estimates. Values and 
errors of the corresponding free estimates give some indication. 

Let us finally emphasize that we have no real knowledge about the analytic 
properties of exponents when d approaches 2. Therefore the be values could be 
affected by systematic effects. 

The article then is organized as follows: in section |] we summarize a few ideas 
about perturbative expansion at fixed d = 3 dimension and ^-expansion. In 
section |3| we briefly recall the Borel summation method based on a conformal 
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mapping of the complex cut plane. Several new variations of the practical imple- 
mentation of the general method are explained. In section ^| we recall the idea of 
the pseudo-epsilon expansion and introduce the exponents' correlation analysis, 
which consists in eliminating the coupling constant between different exponents. 
Section contains a discussion of the numerical results. Finally the new values 
of exponents slightly affect the recently published results |13[ for the equation 
of state of the 3D Ising model, and we present the new determination in section 
10 (as well as a revised version of e-expansion predictions). 



Table 1 

Critical exponents of the O(N) models from d = 3 expansion (present work). 



N 





1 


2 


3 


~ * 

9m 


1.413 ±0.006 


1.411 ±0.004 


1.403 ±0.003 


1.390 ±0.004 


g* 


26.63 ±0.11 


23.64 ±0.07 


21.16 ±0.05 


19.06 ±0.05 


7 


1.1596 ±0.0020 


1.2396 ±0.0013 


1.3169 ±0.0020 


1.3895 ±0.0050 


V 


0.5882 ±0.0011 


0.6304 ±0.0013 


0.6703 ±0.0015 


0.7073 ± 0.0035 


V 


0.0284 ±0.0025 


0.0335 ± 0.0025 


0.0354 ±0.0025 


0.0355 ±0.0025 





0.3024 ± 0.0008 


0.3258 ±0.0014 


0.3470 ±0.0016 


0.3662 ±0.0025 


a 


0.235 ±0.003 


0.109 ±0.004 


-0.011 ±0.004 


-0.122 ±0.010 


LO 


0.812 ±0.016 


0.799 ±0.011 


0.789 ±0.011 


0.782 ±0.0013 


6 = LOU 


0.478 ±0.010 


0.504 ±0.008 


0.529 ±0.009 


0.553 ±0.012 



Table 2 

Critical exponents of the 0(N) models from e-expansion (present work). 



N 





1 


2 


3 


7 (free) 
7 (be) 


1.1575 ±0.0060 
1.1571 ±0.0030 


1.2355 ±0.0050 
1.2380 ±0.0050 


1.3110 ±0.0070 
1.317 


1.3820 ±0.0090 
1.392 


v (free) 
v (be) 


0.5875 ±0.0025 
0.5878 ±0.0011 


0.6290 ±0.0025 
0.6305 ± 0.0025 


0.6680 ± 0.0035 
0.671 


0.7045 ± 0.0055 
0.708 


T] (free) 
r\ (be) 


0.0300 ±0.0050 
0.0315 ±0.0035 


0.0360 ± 0.0050 
0.0365 ±0.0050 


0.0380 ± 0.0050 
0.0370 


0.0375 ±0.0045 
0.0355 


(free) 
(be) 


0.3025 ±0.0025 
0.3032 ±0.0014 


0.3257 ±0.0025 
0.3265 ±0.0015 


0.3465 ± 0.0035 


0.3655 ±0.0035 


to 


0.828 ±0.023 


0.814 ±0.018 


0.802 ±0.018 


0.794 ±0.018 


e 


0.486 ±0.016 


0.512 ±0.013 


0.536 ±0.015 


0.559 ±0.017 



4 



Table 3 

Critical exponents in the 0(4) models from d = 3 and e-expansion (present work). 





d = 3 


e : free, be 


~ * 


1.377 ±0.005 




g* 


17.30 ±0.06 




7 


1.456 ±0.010 


1.448 ±0.015 , 1.460 




0.741 ±0.006 


0.737 ±0.008 , 0.742 


V 


0.0350 ± 0.0045 


0.036 ±0.004, 0.033 





0.3830 ± 0.0045 


0.3820 ± 0.0025 


a 


-0.223 ±0.018 


-0.211 ±0.024 




0.774 ±0.020 


0.795 ± 0.030 


e 


0.574 ±0.020 


0.586 ±0.028 



2 Renormalized </> 4 field theory: e-expansion and 3D perturbation se- 
ries 

In this article the general framework is the (</> 2 ) 2 , O(N) symmetric, quantum 
field theory whose bare action is: 

W) = J {\[d^{x)] 2 + iA 2 ^ 2 (x) + i [ A 4 [^ 2 (x)] 2 }d d x. (2.1) 

We recall that near the critical temperature T c A 2 is a linear measure of the 
temperature. If we denote by A 2c the value for which the theory becomes massless 
(T = T c ) then the parameter t 

t = A 2 - A 2c ex T - T c , (2.2) 

characterizes the deviation from the critical temperature. 

The (</> 2 ) 2 field theory is renormalizable in four dimensions, and to eliminate 
UV divergences (for d < 4 the theory is super-renormalizable) one introduces 
renormalized correlation functions. This involves choosing a renormalization 
scheme and then trading the bare parameters A 2 ,A 4 for a (scheme-dependent) 
renormalized mass m and dimensionless coupling g. The mass parameter m 
is proportional to the physical mass, or inverse correlation length, of the high 
temperature phase. It behaves for t oc T — T c — > 0+ as m oc t v , where v is the 
correlation length exponent (see [f§ for details). 

Renormalization group (RG) arguments tell us that the long distance proper- 
ties of the massless (critical) theory are governed by non-trivial IR fixed points 
g*, solution of the equation 



/3(<7*)=0, with/?V)=u;>0. 
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The anomalous dimensions 77(a) and 172(g), of the renormalized field <p r = 
0/ \[Z and of the renormalized composite operator [</> 2 ] r = (Z^/Z)^) 2 respectively, 
evaluated at g = g* then yield the two independent combinations of critical ex- 
ponents (e.g. r\ = r)(g*)). The explicit forms of the RG functions (3(g), 17(g), 172(g) 
depend on the specific renormalization scheme. 

The space dimension relevant for statistical physics is d = 3 (occasionally 
d = 2). In this case one faces a serious problem: ordinary perturbative expansion 
in g in the massless theory is IR divergent for any fixed dimension d, d < 4. 
A solution to this problem was first provided by Wilson-Fisher's e = 4 — d- 
expansion. The idea is to avoid IR problems by expanding in e = 4 — d as well 
as in the coupling constant g. IR singularities are then only logarithmic and can 
be dealt with. The expansion to the highest order presently available have been 
performed within the minimal subtraction MS scheme, [11, [12 . In this scheme 



the (i-dimensional RG beta function takes the exact form 
The fixed point equation 

^(MS) (SmS' C ) = ' 

can be solved in the form of an ^-expansion. The L-loop expansion of the (3- 
function then yields g±-^ up to order e L . By replacing gr^g by g^ in the per- 
turbative expansion of the anomalous dimensions rj, 172 one finally obtains the 
e-expansion of critical exponents. Note that while g* is scheme-dependent, the 
^-expansion for universal quantities is scheme-independent. 

While this method directly yields a formal expansion for exponents a practical 
problem arises when one wants to determine exponents for a physical value of 
e like s = 1 (d = 3). Indeed the e-expansion is divergent as has been first 
empirically noted in and later confirmed by the large order behaviour analysis. 
A summation method is therefore required to obtain accurate results. 

Following Parisi's suggestion || perturbation series have also been calculated 
directly in three dimensions in the framework of the massive renormalized theory 
where correlation functions of the renormalized field ^ r are fixed by the 
normalization conditions 

T[ 2) (p;m,g)=m 2 +p 2 + 0(p 4 ) , (2.3a) 
( Pl = 0;m,g) = mg. (2.36) 

One may be surprised by the introduction of coupling and field renormalizations 
in a super-renormalizable theory. The reasons are simple, the bare coupling 
constant becomes infinite when the physical mass goes to zero. Simultaneously 
the field renormalization also diverges (see HI]]). 

Series up to six loops obtained in this scheme in Ref. [|7j for N = • • • 3 have 



been generalized in Ref. [|36| to any N. Only recently in [|6| the results for 17 
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and 772 (but not for /3(g)) have been extended to seven loops for N = ■ • ■ 3 (see 
Appendix |A1| ). One problem here is that the value g* of the fixed point coupling 
is affected by summation errors on the /^-function. Errors on g* then induces 
systematic errors for all critical exponents (see section |3|). 

3 Series summation 

Perturbative quantum field theory generates divergent series. Summing such 
series by simply adding successive terms is meaningful only as long as coupling 
constants remain small enough (like in QED). Here however the expansion pa- 
rameter, the fixed point value g* , is a number of order 1: one therefore faces the 
problem of evaluating the sum of divergent series in a non-trivial regime. 

In this article the Borel-Leroy transformation has been used, followed by a 
conformal mapping []33| (a new version of the method developed in || for critical 
exponents) to sum the series. We recall that the Borel summability of the </> 4 
theory in two and three dimensions has been established in [[32| . 

Let S(z) be any (Borel summable) function whose series has to be summed. 
We transform the series: 

S(z)=Y J S k z\ (3.1) 

k=0 



into: 



fc=0 

with: 



POO 

S(z) = J2 B k(b) t b e- t u k (zt)dt, (3.2) 
Jo 



vO + as- 1 
Vl + as + 1 

The coefficients Bk are calculated by expanding in powers of z the r.h.s. of 
equation ( |3.2| ) and identifying with expansion (|3.1| ). The constant a has been 
determined by the large order behaviour analysis. The explicit values are 

a = 0.147774232 x — ^— , (3.4) 
N + 8 v ; 

for the perturbative expansion in d = 3 dimensions and 

3 (3-5) 



N + 8 



for the e = 4 — d expansion. We map the Borel plane, cut at the instanton 
singularity s = — 1/a, onto a circle in the it-plane in such a way to enforce 
maximal analyticity and thus to optimize the rate of convergence (for details see 
e.g. I)- 
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Additional technical details. Following an idea introduced in || for the summa- 
tion of the e-expansion we have in addition made a homographic transformation 
on the coupling constant z to displace possible singularities in the complex z- 
plane: 

z = z'/(l + qz'). (3.6) 

We have looked for values of the parameters q and b for which the results were 
specially insensitive to the order k: in practice the absolute differences of results 
corresponding to three successive orders have been minimized. When several 
solutions were found the less sensitive solution was chosen. Moreover the value 
of b had to stay within a reasonable range around the value predicted by the 
large order behaviour. 

For each series S(z) we have applied the summation procedure both to S and 
1/S. Finally we have introduced "shifts", for each series summing 

S.(z)= (s(z) -J^SkA lz\ 

\ fc=0 / 

In practice only the cases s = (no subtractions) and s = 1, 2 have proven 
useful. Thus for each exponent we have obtained six results whose spread gives 
an indication of summation errors. 

In some examples (in particular g*) shifts have produced strongly oscillating 
results. It has appeared that it would be useful to somehow interpolate between 
shifted series. An idea, new to this work, has been to consider the combination 

S' = (l + rg)S, (3.7) 

where r has been used as a third variational parameter (dividing of course the 
final result by the factor (1 + rg)). The precise value of r has been obtained 
by minimizing the dependence in b. This new additional parameter has proven 
quite useful: it has allowed, as expected, to obtain series with better apparent 
convergence as well as better general consistency. It has also revealed that in 
a few cases the apparent convergence at r = was deceptive, the results being 
unstable with respect to a variation of r. These cases had already been singled 
out by the extreme values of the optimal 6, q parameters. 

The main consequence of this new approach has been a decrease in the values 
of g* though the length of the series has not changed (better agreement between 
different shifts at nonzero r, see Fig.l), and of 7 for N = and d = 3 (best r = 
values of the exponent revealed to be unstable). 

Errors. The summation error for any quantity S(z) has been estimated by 
looking at differences between successive orders, sensitivity to the parameters and 
spread between all results concerning the same exponent (this has also involved 
checking scaling relations). 
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In the case of the 3D perturbative expansion, the total error for each exponent 
S is the sum of the intrinsic summation error at fixed g*, AS, and the error 
induced by the error in g*, Ag*\ 

S = S*±AS±(^] Ag*. (3.8) 
\dg) r 

We thus also give the derivatives of exponents with respect to g*. Two derivatives 
are displayed in Table 4 , all other ones can be deduced from scaling relations. 
The reader can thus infer the sensitivity of exponents to a change in the values 
and errors of g*. 



Table 4 

Critical exponents: Sensitivity to g* determination. 



N 





1 


2 


3 


4 




0.10 


0.18 


0.28 


0.39 


0.50 


dv/dg* 


0.069 


0.11 


0.17 


0.22 


0.29 



Let us stress here that we quote in our tables the total combined error (as well 
in [P]||) while only the intrinsic summation error is reported in the Table 10 for 



the alternative result of 46 



For what concerns the e-expansion the total error is directly given by the 
intrinsic summation error of each exponent and the situation is in principle more 
favourable: the only problem then is that the available series are shorter (they 
are technically more difficult to obtain) and the summation error is then bigger! 

Remarks. The comparison between results coming from direct d = 3 series 
and e-expansion is useful not only to test the accuracy of our numerical methods. 
Their consistency is also important to test various assumptions or properties. 

In the case of the d = 3 expansion we assume more analyticity in the Borel 
plane as has been rigorously proven. Semiclassical instanton analysis indicates 
that our assumption is quite plausible but this is not a proof. Moreover several 
authors (see e.g. |T5|,[7D]]) have argued that RG functions are not regular at 
g = g* . We have of course checked that these singularities, if they exist, are weak. 
Numerical evidence is that all RG functions are at least differentiate at g = g* 
(including (3'(g) which yields u). We cannot of course exclude the situation 
where these singularities are so weak as to escape detection, but strong enough 
to influence results at the level of accuracy at which exponents are determined. 
Our apparent errors could then be underestimated. Nevertheless it should be 
emphasized that if the hypothesis of analyticity in the cut Borel plane holds, 
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the Borel summation should anyway converge asymptotically even in presence 
of confluent singularities. 

For what concerns the e-expansion problems are more serious, since Borel 
summability has not even be proven. Moreover there are indications that UV- 
renormalon singularities could prevent Borel summability, These singular- 



ities are related to the large momentum behaviour of renormalized perturbation 
theory (the "Landau ghost" problem). A plausible conjecture is that quantities 
related to the massless theory are renormalon-free since they can be calculated 
in the theory with UV cut-off. This in particular applies to critical exponents. 
Instead the question remains open for quantities only defined in the massive 
renormalized theory, like the fixed point coupling constant g*(e) defined by (2.3.) 
Note that, because the e series are rather short, empirical evidence is weak. 



4 Pseudo-epsilon expansion and exponents' correlation analysis 

In P|, a method was introduced to try to circumvent the problem of systematic 
errors induced by an error in the determination of g*\ the so-called pseudo- 
epsilon expansion. The idea is to mimic the ^-expansion and introduce a new 
parameter e in terms of which g* is expanded as well as all critical exponents. 
The d = 3 beta-function in the scheme Eqs. ( J2.3[ ) has the form: 

P(9) = -9 + 02{g) 

where @2 begins at order g 2 with a positive coefficient of order 1. We then replace 
the /^-function by a new function (3(g, e) 



and expand g*(e), the solution to 0(g, e) = 0, in powers of e. Eventually we have 
to sum the series for the value of e = 1 to recover the initial equation. 

This method has been systematically used in ||, and this explains why some 
of the new values of exponents we obtain in this work differ less from the previous 
values of ]2[ than the change in g* would lead to expect. 

To apply the same method here, we face the problem that the series for the f3- 
function have not been extended to seven loops, and therefore for the exponents 
7 or u, for example, the information of the additional seven loops term cannot 
be used. However since rj(g) starts only at order g 2 , to determine r\ at loop order 
L, g*(e) is required only at loop order L — 1. This also applies to the exponent 5 
which only depends on 77, and to which we have equally applied the summation 
procedures. 

It follows that for N = 0, 1,2,3 the pseudo-epsilon expansion yields genuine 
seven loop information on 77, 6, with apparent errors much smaller compared to 
six- loop results. 



10 



Still to try to circumvent the problem of shorter /3(g) series and g* determina- 
tion, we have in this work introduced another idea. We have directly eliminated 
the coupling constant between a pair of independent exponents. For example we 
have inverted the relation g \— > 2 — 1/v 

g(v) = Y,9k(2-l/v)\ 

and then expressed other exponents as series in 2 — 1/v. In this way we have 
obtained correlation curves between exponents, which all eventually can be trans- 
lated into relations rj(v). We have applied the same idea starting from the expo- 
nents 7, /?, expanding in powers of 1 — I/7 and 4 — 1//?. 

With this in mind it is interesting to consider the derivatives dr//dv at the fixed 
point, which we thus display in Table 5 . Other derivatives can be deduced, using 
scaling relations. The correlation line can be fixed by taking a point from the 
list of Table 1 . 

Finally let us note that we can push this idea up to expanding the RG (3- 
function in powers of for example 2 — 1/v and solving directly the fixed point 
equation (3(v*) = 0. 

We have tried this idea but the main problem we have faced is that the general 
structure of series generated by this set of transformations is rather complicated 
and therefore the apparent errors are quite large (a problem which already limits 
the accuracy of the pseudo-epsilon expansion). Therefore the method has mainly 
be used as a check of consistency among the data generated by more direct 
summation. It is possible that more accurate constraints could be obtained with 
more work to better understand the convergence of these new series, but we have 
eventually generated so many data that it became difficult to analyze all of them 
with the same care. 



Table 5 

Critical exponents: Correlation between exponents at the fixed point. 



N 





1 


2 


3 


4 


dr]/dv (d = 3) 


0.83 


0.59 


0.43 


0.32 


0.27 



5 Numerical results 

Let us first consider d = 3 results. The values of g* have been obtained by looking 
for the zeros of the summed RG function (3(g). The various methods explained in 
section |3] have been used, shifts 0, 1, 2, generating three set of values for each N, 
depending on three parameters 6, g, r. Quoted errors for g* reflect the apparent 
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convergence with the order k (k < 7), the sensitivity of g* to a variation of the 
parameters 5, q, r around optimal values as well as the spread between different 
shifts. As additional checks we have looked for the zeros of the function u(g)/3(g) 
(a rather arbitrary choice with the weak motivation that the derivative yields 
the exponent 9 = uv) and calculated g* from the pseudo-expansion (see for 
details). Final results are reported in Table 1 (Table 3 for N = 4). 

For what concerns the values of exponents, we have summed the seven (six 
for N = 4) loop series at the values of g* determined before. We have summed 
independently the five exponents 7, /3, u, 5, rj by using the three parameters b, q, r 
for each exponent (and its inverse) and shifts 0, 1 (shift 2 was considered only as 
a check). Again, errors have been estimated by decreasing the order and looking 
at the spread between summation of different equivalent series, as explained 
in section |||. Additional checks have been derived from pseudo-expansion and 
exponents' correlation analysis, introduced in section Table 1 and Table 
3 report the results of the analysis. 

For what concerns the e-expansion the procedure is the same as for the d = 3 
series, apart from the fact that the g* step is bypassed, the series being summed 
at e = 1 for the physical dimension three. 

More precisely we have summed the genuine e series for the exponents (called 
"free" in Table 2 ) and we summed as well the modified series 

in which for each exponent is imposed the exact value at d = 2 (referred to 
"be", with boundary conditions, in Table 2 ). For N = 0, 1 the d = 2 exact 
exponents are obtained from the underlying conformal theories, for N = 2 from 
the identification with the Kosterlitz-Thouless transition, while for N > 2 the 
behaviour near d = 2 can be obtained from the O(N) non-linear cr-model. 

The general conclusions are the following: by imposing boundary conditions, 
we decrease the apparent errors for A = 0. For N = 1 apparent errors remain 
about the same but the central values are slightly modified. For N > 2 the 
convergence of the series with boundary condition is worse and we report in 
Table 2 only central values for the exponents for which the convergence seems 
reasonable. Errors can approximately been inferred from the difference with the 
free values and the corresponding apparent errors. 

6 Updated values for the N = 1 equation of state and critical exponents 



The new values of the critical exponents obtained in this work directly affect the 
determination of the scaling equation of state for the N = 1, d = 3 case, by the 
method presented in [13]. We thus report here the new estimates. The results 
for the e-expansion have been revised too (in particular the errors on amplitude 
ratios have been reconsidered). 
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We recall that our starting point was an estimation of the values of coefficients 
Fk of small magnetization expansion for the derivative of the effective potential 
V (free energy) with respect to the scaled renormalized field z (magnetization) 

^ = z+^ + J2F 2k+1 (g)z^. (6.1) 

k=2 



These coefficients Fk have been summed in |T3| by using the available series up to 
five loops p5| , |26| , p7| and are reported in Table 6 , compared with results of other 



techniques (a misprint in the last digit of the value of g* has been corrected). 
A uniform approximation for the equation of state has then be provided by the 
determination of the auxiliary function h(9) defined by the reparametrization 



(see also [||,||,||]): 

z = P e/{i-e 2 f (6.2) 

h{9) = p- 1 (1 - 9 2 f & F (z(d)) . (6.3) 

The Order Dependent Mapping technique ||67|1 . has been used to improve con- 
vergence of the small ^-expansion by an optimal choice of the parameter p. 
The new result coming from the revised values of 7, f3 is 

h(9)/9 = 1 - 0.762(3) 9 2 + 0.0082(10) 6 4 , (6.4) 

that is obtained from p 2 = 2.86. This expression of h(9) has a zero at 

9l = 1.33 , (6.5) 

to which corresponds the value of the complex root zq of F(z), \zq\ = 2.80 (the 
phase, given by Eq. (|6.2j ), is —in (3). 

The revised e-expansion estimations of Fk (in Table 6 ) and of critical ampli- 
tudes presented in this paper are obtained using the revised 7, (3 of Table 2 and 
the following expression of h(9) (summed at e = 1): 

h(9)/9 = 1 - 0.72(6) 9 2 + 0.0136(20) 9 4 . (6.6) 

It should be emphatized that while critical amplitudes and the equation of state 
are universal quantities, h(9) is not a universal function; in particular the variable 
9 of e-expansion should not be identified to the corresponding variable of the 
d = 3 analysis, because they are defined from a different mapping Eq. (|6~2|), 
(different p and f3). It follows that h(9) of the two methods (and their errors) 
cannot be compared directly. 

Our value of g* from e-expansion (in Table 6 ) has been obtained from our 



own analysis of 0(e ) series of |f70"| . We report in Table 6 also the recent results 
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of fl3~8fl , obtained by a direct summation of 0(e s ) series for Fk (improved by 
imposing boundary conditions at smaller dimensions). 

Widom's scaling function f(x) (with /(— 1) = and /(0) = 1) can easily be 
derived by (numerically) solving the following system: 

f(x)=6- s h(6)/h(l) 

f y/p (6.7) 

x 



el)\e 



Table 6 
Equation of state. 



9 



F 7 x 10 4 



F 9 x 10 5 



e— exp.,this work 



e-exp., |7(J[3§| 
d = 3, this work 
d 



3 m 

HT @ 



HT []37] 
HT p] 



mc mi 



MC @3| 



ERG g4] 



ERG |7] 



23.3 

23.4 ±0.1 
23.64 ±0.07 

23.71 
23.72 ± 1.49 
24.45 ±0.15 
23.69 ± .10 
23.3 ±0.5 

24.5 ± .2 
28.9 

20.72 ±0.01 



0.0177 ±0.0010 
0.01715 ±0.00009 
0.01711 ±0.00007 
.01703 
0.0205 ± 0.0052 
.017974 ± .00015 
.0168 ±0.0012 
0.0227 ±0.0026 
0.027 ±0.002 
0.016 
0.01719 ± 0.00004 



4.8 ±0.6 

4.9 ±0.6 
4.9 ±0.5 

10 



5.4 ±0.7 

23.6 ±4 

4.3 
4.9 ±0.1 



-3.3 ±0.3 
-5.5 ±4 
-7± 5 



-2.3 ±1.1 



-5.2 ±0.3 



From Eq. (|6.7| ) and the revised values of the critical exponents we can calculate 
various critical amplitude ratios that are reported in Tables 7 and 8 and are 
compared with other theoretical and experimental results. The reader can find 
all definitions and more details in [O. (See also [plj for a report on the subject) 



Table 8 
Other amplitude ratios. 





Rq 


i?3 


C 4 /C 4 


HT series |37| 
d = 3, this work 
e-expansion, this work 


0.1275 ±0.0003 

0.12584 ±0.00013 
0.127 ±0.002 


6.4 ±0.2 

6.08 ±0.06 
6.07 ±0.19 


-9.0 ±0.3 
-9.1 ±0.6 
-8.6 ± 1.5 
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Table 7 
Amplitude ratios. 





A+/A- 


c + /c- 


R c 


R x 


e — exp., 


23 


1 


0.524 ±0.010 


4.9 




1.67 


£ — exp., this work 


0.527 ±0.037 


4.73 ±0.16 


0.0569 ±0.0035 


1.648 ± 0.036 


d = 3, m 






0.541 ±0.014 


4.77 ±0.30 


0.0594 ±0.001 


1.7 


<i = 3, this work 


0.537 ±0.019 


4.79 ±0.10 


0.0574 ±0.0020 


1.669 ±0.018 


HT |9,3C 


1 




0.523 ±0.009 


4.95 ±0.15 


0.0581 ±0.0010 


1.75 


MC@ 






0.560 ±0.010 


4.75 ±0.03 






bin. mix. 






0.56 ±0.02 


4.3 ±0.3 


0.050 ±0.015 


1.75 ±0.30 


liqu. — vap. 




0.48-0.53 


4.8-5.2 


0.047 ±0.010 


1.69 ±0.14 


magn. syst. 




0.49-0.54 


4.9 ±0.5 







7 Conclusions 

Before discussing our results, let us review the results for critical exponents 
obtained by other theoretical methods or as well as experiments. 

The previous most accurate determinations of the critical exponents of the 
O(N) vector model, from quantum field theory and renormalization group, have 
been reported in JI],[| and are shown in Table 9 (we refer to these results as 
LG-ZJ). In Table 10 we report the Murray-Nickel (M-N) predictions (direct fit 
of g series) with the authors' preferred choice of g* (they report only summation 
errors; errors from g* should be added). In Table 11 we list some values for 
N = 4 obtained from Pade Borel summation of d = 3 series up to six loops (see 



39] where results for many values of N > 3 are given). An analysis based on 
Order Dependent Mapping |67| of d = 3 series can be found in |]69|] . In Table 
12 we report the previous analysis of e-expansion while in Table 13 we quote some 
recent results. Other available theoretical predictions come from the analysis of 
High Temperature (HT) series in lattice models, Table 14 , and Monte-Carlo 
(MC) simulations Table 15 . Finally in Table 16 we report for completeness 
some estimates from the truncated "Exact Renormalization Group" approach. 



Table 11 

Critical exponents: results of Pade Borel summation for N = 4, J3^ . 



N, Ref . 


g* 


7 


V 


V 


4,1 


36 


1 


1.369 


1.449 


0.738 


0.036 
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Table 9 

Estimates of critical exponents in the O(N) symmetric (0 2 )§ field theory LG-ZJ. 



N 





1 


2 


3 


? 


1.421 ±0.008 


1.416 ±0.005 


1.406 ±0.004 


1.391 ±0.004 


7 


1.1615 ±0.0020 


1.2405 ±0.0015 


1.316 ± 0.0025 


1.386 ± 0.0040 


V 


0.5880 ±0.0015 


0.6300 ±0.0015 


0.6695 ± 0.0020 


0.705 ± 0.0030 


V 


0.027 ±0.004 


0.032 ±0.003 


0.033 ± 0.004 


0.033 ± 0.004 


(3 


0.302 ±0.0015 


0.325 ±0.0015 


0.3455 ± 0.0020 


0.3645 ± 0.0025 


a 


0.236 ± 0.0045 


0.110 ±0.0045 


-0.007 ±0.006 


-0.115 ±0.009 




0.80 ±0.04 


0.79 ±0.03 


0.78 ± 0.025 


0.78 ±0.02 


e 


0.470 ±0.025 


0.498 ±0.020 


0.522 ±0.0018 


0.550 ±0.0016 



Table 10 

Critical exponents: direct fit of d = 3 series (error from g* is not reported). 



N, Ref . 


~9* 


7 


V 


V 


0, [ 

1, [ 

2, [ 

3, [ 


16 
16 
16 
16 


1 


1.39 
1.40 
1.40 
1.39 


1.1569 ±0.0004 
1.2378 ± 0.0006 
1.3178 ±0.0010 
1.3926 ±0.0013 


0.5872 ± 0.0004 
0.6301 ± 0.0005 
0.6715 ±0.0007 
0.7096 ± 0.0008 


0.0297 ±0.0009 
0.0355 ± 0.0009 
0.0377 ±0.0006 
0.0374 ±0.0004 



Table 12 

Estimates of critical exponents: previous e expansion results (LG-ZJ). 



N 





1 


2 


3 


7 


1.157 ±0.003 


1.2390 ±0.0025 


1.315 ±0.007 


1.390 ±0.010 


V 


0.5880 ±0.0015 


0.6310 ±0.0015 


0.671 ±0.005 


0.710 ±0.007 


V 


0.0320 ± 0.0025 


0.0375 ± 0.0025 


0.040 ± 0.003 


0.040 ± 0.003 


P 


0.3035 ± 0.0020 


0.3270 ±0.0015 


0.3485 ± 0.0035 


0.368 ±0.004 


to 


0.82 ±0.04 


0.81 ±0.04 


0.80 ±0.04 


0.79 ± 0.04 



For what concerns experimental determinations of critical exponents a few 
significant results are displayed in Table 17 . 

3D series. In general the new estimates displayed in Table 1 are more accurate 
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Table 13 



Estimates of critical exponents from e-expansion t 76/. 



N 





1 


2 


3 


r(o( £ 4 )) 


1.390 ±0.017 


1.397 ±0.008 


1.413 ±0.013 


1.387 ±0.007 


7 


1.1559 ±0.0010 


1.240 ±0.005 


1.304 ±0.007 


1.372 ±0.006 




0.5882 ±0.0011 


0.631 ±0.003 


0.664 ±0.003 


0.699 ±0.004 



Table 14 

Critical exponents for Ising-like systems: HT series. 



N, Ref . 



7 



a 



ijJV 



o,[ Sa 



15J 
16l 



1,0 

1,® 
1,|5| 
1,|9] 

1, p| 

2, [49| 

2, [51| 

3, [49| 
3,[52| 
4,1! | 



1.1595 ±0.0012 
1.16193 ±0.0001 
1.239 ±0.002 

1.2385 ±0.0025 
1.239 ±0.003 

1.2395 ± 0.0004 
1.239 ±0.003 
1.237 ±0.002 



0.588 ±0.001 
0.588 ±0.001 
0.631 ± 0.003 

0.6305 ±0.0015 
0.631 ±0.004 
0.632 ±0.001 
0.632 ±0.003 

0.630 ±0.0015 



1.2385 


± 





0005 


0.6310 ± 0.0005 


1.237 


± 





004 




1.323 


± 





003 


0.674 ±0.003 


1.323 


± 





015 


0.670 ±0.007 


1.402 


± 





003 


0.714 ±0.002 


1.40 


± 





03 


0.72 ±0.01 


1.474 


± 





004 


0.750 ± 0.003 



0.105 ±0.007 
0.101 ± .004 

0.104 ±0.004 
.108 ± .005 



0.57 ±0.07 
0.54 ±0.05 
0.52 ±0.03 



than the previous LG-ZJ results. They are compatible within errors with the 
previous analysis. A closer inspection shows however some significant changes 
which require discussion. 

The main effect comes from the new (and smaller) values of the fixed point 
coupling constant for iV < 3. The changes are a direct consequence of the new 
techniques we have introduced. In the old calculation LG-ZJ had noticed two 
puzzling features: the optimal values of the parameter b were somewhat large, 
compared to what large order behaviour did suggest. Moreover the three shifts 
0, 1, 2 gave strongly oscillating results. 

Introduction of the new parameter r (see Eq. (|3.7|) ) has shown that the old 
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Table 15 
Critical exponents: MC 



N, Ref. 



7 







0, ||,|4] 

1,0 

1, [56] 
1, [58] 

1, [62] 

2, || 

2 'B 

2, [62] 

3, [63] 

3, [62] 

4, [62] 
4,|1 



1.1575 ±0.0006 



1.324 ± .001 
1.323 ± .002 
1.316 ± .003 
1.3896 ± .0070 
1.396 ± .003 

1.476 ± .002 

1.477 ± .018 



0.5877 ±0.0006 
0.58758 ± 0.00007 
0.631 ±0.001 
0.6289 ± 0.0008 
0.625 ±0.001 
0.6294 ± .0009 
0.664 ± .006 
0.670 ± .002 
0.6721 ± .0013 
0.7036 ± .0023 
0.7128 ± .0014 
0.7525 ± .0010 
0.748 ± .009 



0.3269 ±0.0006 



.362 ± .004 



.3836 ± .0046 



.038 ± .002 

.0025 ± .006 
0374 ± .0014 



.042 ± .002 
0027 ± 0.002 
.041 ± .002 
.038 ± .001 



U - OiO -.007 



.44 
.55 ± .06 



.54 ± .08 
.51 ± .11 



Table 16 

Critical exponents: "Exact Renormalization Group" estimates. 



N, Ref. 


9* 


7 


V 


V 


6 = LOU 


1, [ 

2, [ 

3, [ 

4, [ 
1,[ 
1,[ 
1,[ 


44 
44 
44 
44 
47 
35 
72 


1 
1 

1 


1.726 
1.675 
1.619 
1.566 


1.247 
1.371 
1.474 
1.556 


0.638 

0.700 

0.752 

0.791 
0.618 ±0.014 
0.6262 ± 0.0013 
0.625 ±0.007 


0.045 
0.042 
0.038 
0.034 
.054 

0.030 ±0.005 


0.56 ±0.07 
0.48 ±0.04 



apparent convergence corresponded to an unstable region of parameters. By 
varying r we find a region where these problems are solved to a large extent: 
the results are less sensitive, various shifts agree, and all parameters have more 
reasonable values. Figure 1 exemplifies this situation for the N = case. 

Another example exhibited a similar instability atr = 0: 7, iV = <i = 3. 
This (as well as the decrease induced from that of g*) explains the new different 
value we obtained here. 

Finally N = 3 values show a consistent effect: the three exponents 7, f3 
increase. This simply suggests that N = 3 errors had been underestimated in 
LG-ZJ analysis. 
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Table 17 

Critical exponents: selected recent experiments 



N, Ref . 


7 




P 


a 


6 = U3V 


0, 
1, 

1,[74, 

1, 
2,H 

2, 
3, 


i 

p 

78 
79 


79] 


1.25 ±0.01 
1.233 ±0.010 


0.586 ± 0.004 

0.64 ±0.01 

0.6708 ± 0.0004 
0.6705 ± 0.0006 


0.327 ±0.002 


0.107 ±0.006 
0.109 ±0.006 

-.01285 ±0.00038 


0.57 ± .09 
0.51 ±0.03 

0.61 ±0.06 



gO vs r, N=0 

gO 



♦ shift=0 

* shift=l 
- - *- ■ shift=2 



J)* .02 5"" 0.0 5 0.075 0.1 0.125 0.15 

Fig. 1 Values of g* , N = 0, as a function of the parameter r for shifts 0,1,2 

Conversely the values of v, N = 0,1,2 are quite stable, even though the 
corresponding values of g* have changed. The reasons are a very good apparent 
convergence of the pseudo-epsilon expansion (on which previous analysis partially 
relied) at previous order. 

The apparent errors have generally been reduced, as should be expected, ex- 
cept for N = 3 (see the comment above) . The improvement is specially significant 
for the exponent r\ that was poorly determined before. With a few exceptions the 
general trend for a given exponent is the increase of the direct summation error 
with N. This effect has a simple explanation: our summation method relies on 
the large order behaviour analysis, and the asymptotic regime sets in later when 
N increases [^0| . Perhaps a clever use of the knowlege coming from the large N 
expansion could improve the situation. 
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Note finally that the term added to two of the three RG series has allowed not 
only to decrease apparent errors but also to estimate them more reliably. 

For several exponents (such as e.g. 7(A = 1)) errors are now dominated by 
errors induced by the determination of g*. To further improve the situation it 
will be necessary to also add a new term to the RG /3-function. 

(Free) e-expansion and 3D series. Comparison with the previous LG-ZJ e esti- 
mates shows no striking effect, small deviations being due to the use of corrected 
series. 

For the exponent v the consistency between 3D and e results remains very 
good at all N. The situation has markedly improved for the exponent 7, N = 0: 
there is still a systematic discrepancy in the central values for the exponents 
(about 0.002), but the difference is reduced by more than a factor two, which 
is quite encouraging. A similar comment applies to the central values of the 
exponent r\ N = 0, • • - , 3, where the discrepancy is also reduced by a factor two. 
For what concerns the N = 4 prediction, the agreement with the corresponding 
d = 3 results is quite satisfactory but apparent errors are large. 

One point should however be stressed: since the series are shorter it is more 
difficult to assess apparent errors and the errors we quote are thus less reliable 
than for the d = 3 series. 

Free and be e-expansion. For the £-expansion we report a second set of values, 
obtained by imposing the exact d = 2 values, referred as be (i.e. with bound- 
ary conditions) in Table 2 to distinguish them from the unconstrainted values 
denoted by free. Some remarks are in order about the be approach. First we 
do not know the analytic structure of exponents near d = 2. The only piece of 
evidence comes from using the e-expansion for d = 2. For N = the agreement 
with exact results is quite good [|], suggesting that d = 2 is a regular point. For 
N = 1 the agreement is less striking, leaving room for some complex behaviour. 
For values N > 2 there are indications (IR-renormalons of the d = 2 non-linear 
(T-model) that d = 2 corresponds to an essential singularity (probably non-Borel 
summable). It is then easy to construct examples where fitting leads to worse 
results. 

Second in the case N > 2 the new series have a more complicated structure 
which makes summation (which is a form of extrapolating series to higher orders) 
and even more error estimation difficult. 

Although the two analyses give compatible results, it happens that be has the 
general tendency to give values of 7, v for N > 2 larger than free ones and values 
of 77 for N > 2 smaller than free ones. This point is understood by remembering 
that at d = 2 the corresponding 7,^,1/77 are infinite and thus it is reasonable 
that imposing this behaviour at d = 2 tends to increase the value of the exponent 
at d = 3. 

Finally it is also remarkable that for A = the be errors are smaller than the 
free ones. 
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Other methods. The general agreement between the HT series, the MC results, 
and the new d = 3 determinations is in general improved. This in particular ap- 
plies to the SAW where recent long simulations provide very accurate estimates. 



NU N = 2 



O . 6 7 6 

O . 6 7 4 

O . 672 

O . 6 7 

O . 6 6 8 

O . 6 6 6 

O . 6 6 4 



" T T 



♦ 


He 


1 ± 




He 


go 




ci3 


gz 




ci3 


1 z 




ep 


gz 


♦ 


ep 


1 z 


* 


ep 


pv 




mc 


□ 




mo 


to 




ht 


to c 



Fig. 2 Comparison between various estimates of the exponent v, N=2: He li from 
[14, 77 1, He go from J78| , d3 gz from present work (d = 3), d3 Iz from LG-ZJ (d = 3), 
ep gz from present work (e), ep Iz from LG-ZJ (e), ep from j7(| (e), mc j from fl6l| , 
mc b from [62], ht be from [p9|. 



Experiments. The improved agreement of present results with the recent mea- 
sures on superfluid Helium systems, iV = 2, is remarkable and is displayed graph- 
ically in Figure 2. In spite of our efforts the best experimental value is still much 
more accurate than the theoretical estimate. In all other cases the agreement 
with experiments is good as seen from Table 17 . 
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APPENDICES 



Al Series d = 3 

We report here the seven loop series for the O(N) symmetric (4> 2 ) 2 theory, N = 
0,1,2,3, computed by Murray and Nickel [||§. The functions rj(g) 112(g) below 
are defined by 

dlogZ d\ogZ 2 

v(g) = m — A — V2(g) = m — A — 

dm dm 

where the renormahzation constant are defined by <po = y/Z(f> r and (4> 2 ) r = 
^(4>o) 2 (subscript 0,r indicate respectively bare adn renormahzed fields). The 
critical exponents r], v can be found by the identification 77 = r\(g*) and v = 
(2 + 112(9*) — ??(<?*)) -1 - The symmetry number N is reported in square brackets 
below. 



a 2 

77 [0] = f— + 0.0007713750<7 3 + 0.0015898706^ 4 - 0.0006606149^ 5 
108 

+ 0.0014103421# 6 - 0.001901867^ 



80 2 

r)[l] = — — — h 0.0009142223# 3 + 0.0017962229<? 4 - 0.0006536980s 5 
729 

+ 0.0013878101<7 6 - 0.001697694^ 



80 2 

v [2] = — + 0.0009873600# 3 + 0.0018368107^ 4 - 0.0005863264# 5 
675 

+ 0.0012513930^ 6 - 0.001395129^ 



r)[3] = ^jL + 0.0010200000<? 3 + 0.0017919257# 4 - 0.0005040977# 5 
+ 0.0010883237^ 6 - 0.001111499^ 7 



~2 

rjzlO] = Z£ + !L - 0.0357672729^ 3 + 0.0343748465# 4 - 0.0408958349# 5 
4 16 



+ 0.0597050472(f - 0.09928487# 7 
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772 [1] = -JL + - 0.0443102531/ + 0.0395195688/ - 0.0444003474/ 
+ 0.0603634414/ - 0.09324948/ 

V2 [2] = + ^- - 0.0495134446/ + 0.0407881055/ - 0.0437619509/ 
+ 0.0555575703/ - 0.08041336/ 

77 2 [3] = d*£. + 1^- - 0.0525519564/ + 0.0399640005/ - 0.0413219917/ 
11 121 

+ 0.0490929344/ - 0.06708630/ 



